On the Theory of Space Charge Between Parallel 
Plane Electrodes 

By C. E. FAY, A. L. SAMUEL and W. SHOCKLEY 

The problem of the potential distribution, current, and electron 
transit time resulting from the perpendicular injection of electrons 
into the space between parallel planes is considered. The elec- 
trons are assumed to be injected uniformly with velocities cor- 
responding to the potential of the plane through which they are 
injected. Consideration of all possible solutions of the basic equa- 
tion shows that four general types of potential distribution are 
possible. Curves are given which enable the easy calculation of 
transmitted current and transit time and show the complete po- 
tential distribution for any concrete example. The case for cur- 
rent injected through both planes is also considered. 

The complete mathematical treatment is given in the appendix. 

CPACE charge has been studied extensively since the publication of 
^ the first papers on the subject by Child 1 in 1911 and Langmuir 2 
in 1913. These papers in common with many which followed 3 - 4 - 5 - 6 - 7 
dealt for the most part with potential distributions which occur when 
electrons are injected into a region with relatively small initial velocities. 
The problem of space charge between parallel planes when the 
electrons possess arbitrary initial velocities, although contained 
implicitly in some of this early work, was first considered in detail by 
Gill 8 in 1925. Gill appears to have clearly understood the phe- 
nomenon but did not publish a complete analysis. Other workers, 
beginning with Tonks 9 in 1927, have considered various aspects of 
the problem, 10, u and recently Plato, Kleen, and Rothe 12, 13 published 
an extensive analysis. They did not include transit time calculations 
and their published curves are not easily adaptable for numerical 

I C. D. Child, Phys Rev., 32, 492 (1911). 

■ I. Langmuir, Phys. Rev., 2, 450-486 (1913); also Phys. Zeits, 15, 348 (1914). 

3 W. Schottky, Phys. Zeits, 15, 526 and 624 (1914). 

I P. S. Epstein, Vehr. d.D. Phys. Ges., 21, 85 (1919). 

8 T. C. Fry, Phys. Rev., 17, 441 (1921); Phys. Rev., 22, 445 (1923). 
I. Langmuir, Phys. Rev., 21, 419 (1923). 

7 1. Langmuir and K. B. Blodgett, Phys. Rev., 22, 347 (1923); Phys. Rev., 24, 49 
(1924). 

8 E. W. B. Gill, Phil. Mag., 49, 993 (1925); Phil. Mag., 10, 134 (1930). 

9 L. Tonks, Phys. Rev., 30, 501 (1927). 

10 H. C. Calpine, Wireless Engineer, 13, 473 (1936). 

II Myers, Hartree and Porter, Proc. Roy. Soc. 158, 23 (Jan. 1937). 

»G. Plato, W. Keen, and H. Rothe, Zeit.f. Phys., 101, 509 (July 1936). 
18 H. Rothe and W. Keen, Telefunken Rohre No. 9 (April 1937). 
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calculations. In this country Salzberg and Haeff 14 have presented a 
solution, apparently quite complete, which has not as yet been pub- 
lished. The present independently derived solution differs from other 
published work in the manner of presentation and in an attempt at 
completeness. The results are presented in the form of curves which 
give potential distributions, electron transit times, and current-voltage 
relations, directly in terms of units which have a simple physical 
significance. 

Statement of the Problem 

In common with some of the earlier treatments, the present dis- 
cussion will be restricted to a consideration of the steady state potential 
distributions which can exist in an evacuated region between two 
parallel planes at known potentials when electrons having normal 
velocities corresponding to these potentials are injected into the region 
through one or both planes. For definiteness the case in which all 
of the current is injected through one plane will be considered first. 
It will then be shown that the results may be applied to the more 
general problem. It is supposed that the electrostatic potential is a 
function of one rectilinear coordinate only and that all the electrons 
move parallel to this coordinate and have the same total energy. 
The assumption of equal energy electrons leads to indeterminancy at 
planes where the potential and its gradient are both zero. At these 
planes the existence of non-uniform electron velocities will be recog- 
nized in so far as it provides a selective mechanism to resolve the 
mathematical indeterminancy. 16 

Units 

It has been found convenient to express voltages, currents, distances 
and transit times in terms of some derived units which are related to 
these quantities in certain ways by the simple Child's law equation 
for space charge limited current and by the corresponding transit 
time equation. The use of such derived units makes it possible to 
present a limited number of curves which are then applicable to a 
wide variety of conditions and may be used with a minimum amount 
of computations. The physical significance of the units also simplifies 
the interpretation of the results. 

In order to understand how these units are obtained consider the 
hypothetical situation shown in Fig. 1 in which an electron current 
from a space charge limited cathode is injected into the region to the 

11 B. Salzberg and A. V. Haeff, Proc. I.R.E., 25, 546 (May 1937).— Abstract only. 
This paper appears in full in the January 1938 issue of the R.C.A. Review. 

« This procedure has been justified by Langmuir. See Reviews of Modern Physics, 
vol. 3, pp. 237-244 (1931). 



SPACE CHARGE BETWEEN PARALLEL PLANE ELECTRODES 51 



right of a plane at a potential V\. u From Child's equation the current 
in amperes per square centimeter is given by 

7,3/2 a 2 y 3/2 

/ = 2.33 X 10- fl -^- = 2_I*_ ( am P s - P er SQ- cm -) (1) 



or solving for So 



V,3/4 n J/,3/4 

5 = 1.527 X 10- 8 -^r = ^- (centimeters). 



(2) 



Jl/2 Jl/2 

Accordingly, whenever the conditions at the first plane are represented 
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Fig. 1 — Hypothetical conditions to assist in visualizing the unit of distance So. 



by J and Vi, distance from the first plane may be measured in units 
of S . The distance in centimeters is then 



5 = <rSo. 



(3) 



Similarly a natural unit of potential is V\ and the potential of any 
plane at a distance S is then given by 



V = tpV\. 
Potential Distributions 



(4) 



All possible potential distributions to the right of the first plane 
can now be expressed in terms of <p as a function of a. The mathe- 
matical derivation is straightforward and is given in the appendix. 

"This analogy is useful in getting a physical picture of the units but it must 
not be carried too far as will be evident when the problem of reflected currents is 
considered. 
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The results are shown in Figs. 2, 3 and 4. Figures 6 to 11 refer to 
current voltage relationships discussed in a later section. 

To find the potential distribution when a second plane at a distance 
a away from the first plane is held at a potential <p, one enters the 
figures with the values of <p and a. Any member of the family of 
curves drawn in heavy solid lines passing through this point represents 
a possible potential distribution. The regions occupied by the heavy 
curves are bounded by certain limiting curves identified by lower case 
letters. These curves correspond to certain space charge conditions 
between the planes and are indicated by the same letter whenever 
met. The dashed curves refer to transit time and are explained below. 
Situations corresponding to low potentials and large spacings occur, 
for the most part, on Fig. 2, while those for high potentials and small 
spacings occur on Fig. 3. An overlap region exists for which potential 
distributions may be found on both families of curves. The interpre- 
tation of this is found by inquiring more closely into the nature of the 
curves of Figs. 2, 3 and 4. 

There are in all four distinct types of potential distributions. 
These are : 
Type A — The second plane is at a negative potential. All the injected 

current is reflected. Potential distributions are the same as 

for the case of temperature limited emission with the current 

27 from the zero potential plane. 
Type B — Both planes positive with potential zero between them. 

Injected current partially transmitted and partially reflected. 

Potential distributions corresponding to space charge limited 

emission on both sides of a virtual cathode at the zero of 

potential. 
Type C — Both planes positive with a potential minimum (at a positive 

potential) between them. Complete transmission of injected 

current. 
Type D — Both planes positive with no minimum between them. 

Complete transmission of injected current. 
The curves on Fig. 2 are seen to fall into two groups, those which 
extend to negative values of <p (marked with values of /3) and those 
which contain the value <p = but which remain positive. The first 
group (Type A) obviously corresponds to conditions under which all 
of the injected current returns to the first plane. The slopes of these 
potential curves at the reflection planes are not zero but are related to 
the parameter /3 (shown on the curves) by the equation 

^=-tV2>". (5) 
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Fig. 2 — Potential distributions of the A and B types. The solid lines are the potential curves; 
the broken lines indicate the transit. times. 
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For the second group of curves (Type B) a certain fraction of the 
current denoted by Z (and so indicated on the curves) will be trans- 
mitted while the fraction 1 — Z will be reflected toward the first 
plane. These potential curves have zero slopes at the reflection planes 
so that they correspond to solutions of Child's equation on both sides 
of a so-called virtual cathode existing at the reflection plane. For 
each value of Z from to 1 a possible potential distribution is obtained. 

It will be observed that the portions of the solid curves to the left 
of the zero points are drawn lighter than the rest. These portions 
correspond to potential distributions resulting from conditions with 
reflected current and so while applying as extensions of the heavy 
curves to the right they cannot be entered directly with values of <p and 
<r. Should values used to enter the figure fall in this region the ab- 
sence of a potential zero is indicated and the correct distribution is 
obtained by entering Fig. 3. Before leaving Fig. 2 the existence should 
be noted of a small domain lying between the curves marked b and c 
for which two different B solutions are possible for the same injected 
current, spacing, and potential, one corresponding to larger values of Z 
than the other. The significance of this region will be apparent when 
the current voltage relationships are considered. All solutions of the 
A and B types are represented on Fig. 2. 

Solutions of the C and D types are to be found on Fig. 3 except for a 
small overlap region which for clearness is shown on Fig. 4. Values of 
<p and a which cannot be entered on Fig. 2 as well as values common to 
both B and C types are to be found on these figures. The additional 
overlap solutions of the C type are shown in Fig. 4. These are ex- 
tensions of the potential distribution curves of Fig. 3 after they reach 
the right-hand boundary curve where they turn inward as shown and 
overlap the other solutions. 

For convenience each curve of the C type is labeled by the value of 
its minimum potential. The parameters of the D curves do not have 
this simple physical significance. However, for all cases the value of 
the curve parameter is simply related to the electric field at the first 
plane by the equations: 

typeC ^ = -IVl-^.-n. 1 ' 2 . (6) 



type D 



d f a = |vn^, (?) 

d<p 



7=±|VT+^. (8) 

da 5 

dip 4V2 



type A ^ = - 1 T :Vr + ^- (9) 
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3 — Potential distributions of the C and D types. The solid lines are potential curves, the 
broken lines indicate the transit times. 
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CT = DISTANCE IN UNITS OF S (So= 1.53 X I0 -3 ^Tr) 

Fig. 4 — Potential distributions of the C overlap type. 
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It should be noted that while some of the curves in the D region 
shown in Fig. 3 resemble the portions of the curves drawn in light 
lines on Fig. 2, the curves in Fig. 3 have essentially different physical 
content, representing solutions in which all of the injected current is 
transmitted. Values of a and <p may, therefore, be entered in all 
regions on Fig. 3 while entering values on Fig. 2 are restricted to the 
regions covered by the heavy lines. 

The occurrence of overlap regions indicates that for certain boundary 
conditions more than one type of potential distribution may exist. 
In the practical case the choice between these various possible solu- 
tions depends upon the manner in which the boundary conditions are 
established. 17 

Current Injection from Both Sides 

The present analysis although derived on the assumption of current 
injection from one side only is equally applicable to the situation 
shown schematically in Fig. 5. The potential distribution curves 
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Fig. 5 — Schematic representation of the general conditions under which the potential 
distribution analysis may be applied. 

which occur when the currents I a and 2& are injected from opposite 
sides will be identical with those obtained for an assumed injection 
from one side only of a current / equal to the numerical sum of I a and 
lb- Values of Vi, <p and a are chosen to correspond with an assumed 
direction of injection and these values are entered on the figures. 
For solutions of the B type, different distributions will result depending 

"This matter will be treated in more detail in the section on circuit charac- 
teristics. 
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upon the assumed direction of I. In some instances all solutions may 
correspond to physically realizable distributions for the double in- 
jection case. The solutions will, in general, require that the potential 
zero be located at different places and may require that the net flow 
of current across the zero plane be in opposite directions. The 
physical reality of each of these distributions must be checked by 
noting the value of Z for the indicated solution and solving for the 
derived current (2 — Z)I on the injected side of the zero potential 
plane and the derived transmitted current ZI. If these are respec- 
tively greater than the numerical values of I a and h (paired consistent 
with the assumed direction of injection) the indicated solution is 
possible. In addition there may exist a solution in which the currents 
I a and lb are each totally reflected at two different zero potential 
planes, separated by a region of zero potential. The possible existence 
of this solution must be tested for separately. 

Electron Transit Time 

So far no use has been made of the dashed curves shown in Figs. 2, 
3 and 4. These give the electron transit time from the first plane 
to any desired plane. The significance of the unit of time is to be 
found by again referring to Fig. I. The time an electron takes to 
travel from the cathode to the first plane is given by 

T7U/4 V\ 1/4 

/o = c -^ = 7.72 X 10- 11 ^- (seconds). (10) 

This value of to is a natural unit of time to be used whenever the condi- 
tions at the first plane are expressed by Vi and I just as S is a natural 
unit, of distance. Accordingly, electron transit times from this plane 
to any other plane are expressed in units of t - The time in seconds is 
then 

T = T/ . (11) 

Transit times for reflected electrons are not given but may be com- 
puted by taking the time to the reflection plane and adding to it the 
returning time. This latter is obtained by taking the difference be- 
tween the time to the reflection plane and the direct time to the plane 
under consideration. 

Circuit Characteristics 

In many practical cases the first plane will coincide with a mesh or 
grid electrode in a vacuum tube and the second plane will correspond 
with a plate electrode. Under these conditions the current-voltage 
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relations which exist at the second plane are of considerable interest. 
Since the distance between the planes is now fixed, the parameter <r is 
not appropriate. However, the spacing d between planes and the 
voltage Vi serve to define a natural unit of current density to, which 
is indicated by the hypothetical situation of Fig. 6. The value of io is 
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Fig. 6 — Hypothetical conditions to illustrate the significance of the unit of current ** . 



to = a 



T/J/2 



Vl W 



' "-= 2.33 X 10~ 6 ~- (amperes per sq. cm.) 



d 2 



(12) 



The current in units of i is denoted by 7 so that the injected current / 
in amperes per sq. cm. is given by 



I = yi . 



(13) 



Similarly, the transmitted current under conditions corresponding to 
type B potential distributions will be 



ZI = Zyi . 



(14) 



Specifying 7 and <p is equivalent to specifying a and ip. 1 * For this 
reason many of the limiting curves of the circuit characteristic plots 
are simply related to those shown in the potential distribution plots 
and are designated by the same letters. The regions for which the 
different types of potential distributions may exist are shown in Fig. 7 
now in terms of 7 and <p. In Fig. 8 the regions are defined in terms of 
ip and the transmitted current Zy. The boundaries of the C and D 

18 This is a consequence of the relationship y = a 2 , which is derived in the ap- 
pendix. 
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r = INJECTED CURRENT IN UNITS OF 
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regions are of course the same in the two plots (i.e., Z = 1) while the 
boundaries of the B region are changed, since only the portion Zy of 
the injected current is transmitted. The A region, which corresponds 
to negative values of <p for all values of 7 and gives Zy = 0, is not 
indicated. It may be of interest to note that for any given value of <p 
it is impossible to transmit more than a certain maximum current, 
shown by curve a, and that this maximum occurs for a C type distri- 
bution. For injected currents greater than this amount only type B 
distributions are possible. 

An enlargment of the B region is shown in Fig. 9a and a portion of the 
region to a still larger scale in Fig. 9b. Lines for constant injected 
current 7 are plotted with the potential <p and the transmitted current 
Zy as the coordinates. These plots correspond to plots frequently 
used to describe vacuum tube characteristics and may be used ac- 
cordingly. For values of 7 less than four and values of <p less than 
unity, double values of Zy appear. In the region lying between the lines 
b and c on Fig. 9, the slopes of the constant 7 lines are negative corre- 
sponding to conditions which are unstable unless sufficient resistance is 
included in the external circuit. Conditions outside of this region are 
always stable and need no further comment. 

Solutions of the C type are shown in Fig. 10a. For these solutions 
the injected and transmitted currents are equal. The depth of the 
potential minimum between the two planes is, however, of interest 
and is shown in the figure in terms of the lesser of the two boundary 
potentials. For values of <p greater than unity the value of <p min. 
(i.e., the minimum potential in units of V\) is indicated while for 
values of <p less than unity the value of <p m i». (i.e., the minimum po- 
tential in units of V 2 ) is shown. In the enlargement, Fig. 106, both 
values are given. For values of ip m in. = 1, curve e, and <p m in. = 1, 
curve d, the potential minimum is equal to the potential at one of the 
planes and is located at this plane. For currents greater than the 
value indicated by these limiting curves the potential minimum be- 
comes deeper and moves away from the plane while for lesser currents 
it disappears entirely. The minimum value to which the potential 
may be forced by increasing the injected current before the distribu- 
tion changes precipitously to one of the B type is indicated by the 
values on curve a. 

Again to avoid confusion the overlap type C region is shown sepa- 
rately in Fig. 11. Conditions are somewhat different for this region 
in that with any assumed value of <p a minimum potential of (but 
with complete transmission of injected current) occurs for a certain 
value of 7. With increasing injected current, the potential minimum 
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rises, finally reaching the limiting values indicated on the upper curve 
of Fig. 11; the limiting conditions there are precisely the same as for 
curve a of Fig. 10. Further increase of injected current produces a 
transition to type B solutions. 

It seems questionable that this overlap type C condition can exist 
in a practical case. An investigation of this question involving ex- 
ternal circuit considerations is beyond the scope of this paper. 

Transition Between Distribution Types 

The physical choice between the different possible potential dis- 
tributions which may exist with a given set of boundary conditions is 
determined by the sequence in which the boundary conditions are 
established. Extreme values of any parameter are seen from Figs. 7 
and 8 to lie in regions for which only one solution is possible. If the 
boundary conditions are varied slowly and continuously from these 
values, the indicated type of distribution will persist until the limit of 
this region is reached at which time a sudden transition must occur 
to another indicated type of distribution. Inspection of Figs. 7 and 8 
will show that at such transitions only one other type of distribution 
is ever possible. The determination of the correct physical distribu- 
tion can thus be made without ambiguity. 

Certain peculiarities are, however, to be noted. A survey of all 
possible transitions in which y and v> are treated as independent 
variables will indicate that, starting from extreme conditions and 
changing conditions continuously in the same direction, distributions 
of the overlap C type shown in Figs. 4 and 11 never occur. A second 
peculiarity has to do with the unstable region of type B solutions 
shown on Fig. 9. When this region is entered with insufficient 
resistance in the external circuit, instability results with a sudden 
transition to a corresponding stable type of distribution. 

The space model shown in Fig. 12 has been found to be of value 
in visualizing problems involving transitions. The three coordinates 
used in its construction are the second electrode potential <p (to the 
right) the injected current 7 (to the left) and the transmitted current 
Zy (vertical) . Solutions corresponding to potential distributions of the 
C and D type, for which Zy = y, appear as a celluloid plane inclined at 
45 degrees, on which the values of the potential minima are indicated. 
Solutions of the C overlap type have been omitted. Solutions of the B 
type are represented by the concave surface of the model. Viewing the 
model from a different angle as shown in Fig. 13 (where the celluloid 
sheet is removed) the unstable B region appears as an over-hanging 
cliff extending between 7 values of 0.5 and 4. Bounding curve c 
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separating the unstable region from the stable region is drawn on the 
model and lies along the vertical surface of the cliff. The small 
region lying between 7 values of and 0.5, for which direct transitions 
between A solutions and D solutions occur, can be seen in both views 
as a small vertical cliff. 

Solutions corresponding to varying conditions will be given by the 
position of a movable particle confined to the surface of the model. 




Fig. 12 — A three-dimensional model illustrating the current- voltage relationships. 
Solutions of types" C and D appear on the celluloid surface and solutions of the B 
type appear on the concave surface. 

When this particle lies on the celluloid surface it may be moved around 
at will and will not "fall off" unless an attempt is made to go to 
larger values of the injected current 7 than is permitted by the bound- 
curve a (the curved edge of the celluloid). The particle may, however, 
"fall off" at this edge and land on the concave surface from which it 
can escape only by climbing up the surface to the bounding curve b 
which is common to both surfaces. The overhanging cliff is, however, 
particularly treacherous, for unless the particle is supported by the 
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existence of sufficient resistance in the circuit connected to the second 
electrode, it will go through the model coming to rest on the celluloid 
surface above. With no resistance in the external circuit this curious 
behavior will occur at the vertical part of the cliff indicated by the 
bounding curve c. The presence of some resistance will cause this to 
occur for larger values of Zy. The exact place can of course be 
determined by noting the point of tangency between a resistance line 
drawn on Fig. 9 and one of the characteristic curves — just as the 
stable operating conditions for any negative resistance device is found. 




Fig. 13 — A different view of the concave surface of the model shown in Fig. 12. 
(The celluloid surface has been removed.) 

The entire surface is stable if the external resistance is greater than 
approximately 0.25 in units of <p/Zy. 

Illustrative Examples 

/. Class B Operation of a "Critical Distance'' Tetrode 

Since the purpose of "critical distance" operation is to prevent the 
passage of secondary electrons from the plate (our second plane) to 
the screen grid (our first plane) , the existence of a potential minimum 
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is desired ; however, type B solutions with some reflected electrons 
would be objectionable. The permissible operating range must 
therefore lie in the C region on Fig. 7. If, for example, we operate 
the tube with V2 = V\ as the quiescent point and with a resistance 
load, the operating curve is a straight line through the point 7 = 0, 
<p = 1. The upper excursion of this line is set by the upper limit, 
curve a, of the C region. For characteristics which cross curve d and 
enter the D region a loss of suppression would result. Maximum 
efficiency occurs with the largest possible excursion of F 2 , that is for 
a line tangent to the d curve at the quiescent point. The minimum 
instantaneous value of V2 for the conditions specified is seen to be 
25 per cent of its quiescent value as contrasted with a usual value of 
about 10 per cent for a pentode when similarly operated. For the 
instantaneous value of F 2 to fall to 10 per cent of its initial value 
(without loss of suppression) the quiescent point must be at 7 = 0, 
(p = 1.3 with the operating line tangent to curve d at about (p = 0.5. 
Operation in the D region above <p = 1 does not 'result in a secondary 
current from plate to grid. 

77. Triode Operation in the Positive Grid Region 
Figure 1 may be taken as representing an idealized triode. The 
second plane represents the plate so that V2 = V p and the first plane 
represents the grid plane with effective potential V\. We let 

v, = v a + v P i», 

where /z will be a variable if space charge is formed in the plate-grid 
region and may differ considerably from its "cut-off" value. 

If the fraction of the electrons stopped by the grid is negligible and 
if operation does not enter the B region, with disturbing effects 
produced by reflected electrons, then Fig. 1 becomes an accurate 
representation and in view of the choice of units a is a geometrical 
constant of the tube. Hence the operating characteristic is a constant 
7 line (7 = <7 2 ) on Fig. 7. 19 

For normal operation the characteristic must not extend into the B 
region. If the grid is to be driven far positive it is desirable to stay 
in the C region for all values of <p < 1 to prevent secondary electrons 
from leaving the plate. The minimum value of ^ as a function of 7 
may be read from the intersections of the constant 7 lines with limiting 
curve a if 7 > 2 and limiting curve d if 7 < 2. Now since 



<P = 



V P V t 



Fl V. + & 



19 Proof of the relationship 7 = a 1 is contained in the appendix. 
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and 

Vg Vp = = > 

if n and a are known, the maximum value of VJ Vp may be calculated. 
Fig. 7 shows that a value of 7 just greater than 2 will give the lowest 
minimum value for <p (<p = 0.07) by intersection with curve a. In 
this region <p/fi <C 1 so that 

-~- max. = — = 14. 

Vp J <p 

For values of 7 only slightly less than 2 the minimum value of *p rises 
abruptly, being given by the intersection with curve d, while for values 
greater than 2 the rise in <p is not so abrupt. 

It should be noted that for the usual cylindrical triode the departure 
from the parallel plane case is not very marked in the grid-plate 
region while it can "be taken into account by the usual cylindrical 
formula for the cathode grid region. It therefore becomes possible 
to extend this analysis in an approximate sort of way to cover the 
cylindrical case. 

APPENDIX 

Definition of Symbols 

Vi = Potential of the first plane in volts, 
V = Potential of a second plane in volts, 
<p = Ratio of the potential at the second plane to the potential of the 

first plane, 
I = Injected current in amperes per sq. cm., 
x = Distance in centimeters, 
d = Spacing between the two planes in centimeters (used when this 

distance is regarded as a fixed quantity), 
5 = Distance between the two planes in centimeters (used when this 

distance is regarded as a variable), 
a = Distance between the planes expressed in units So, 

So = 1.527 X 10- 3 ^- (centimeters), 

a = 1.527 X 10- 3 = [10 ,3 (2e/wc) 1 / 2 /9 7 rc 2 ] , / 2 , 2 ° 
a 2 = 2.33 X 10- 6 , 
Z = Fraction of the injected current I which reaches the second 

plane, 
/ = Time in seconds, 
T = Electron transit time between the planes in seconds, 
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t = Electron transit time expressed in units of to, 

h = 7.72 X lO" 11 ^- 4 (seconds), 

c = 7.72 X 10-" = 3a(tnc/2eyi*10r*« 

y = Injected current expressed in units of i , 

Vi 3 ' 2 a 2 Vi 3 ' 2 
io = 2.33 X 10~ 6 , 2 = — -ij — (amperes per sq. cm.), 

a = Ratio of the potential at a potential minimum when one exists 
to the potential of the first plane, otherwise merely a con- 
venient constant of integration. (In the text the symbol 
tPmin. is used in the cases where the physical significance can 
be attached.) 

/3 = An integration constant similar to a but associated with an 
opposite sign. 

Relationship between y and <r 
By definition 

_aV l 3 '*_S 



and 

Solving (1) for <r 

Solving (3) for y 



u -—dT-y' (2) 



en/2 

(3) 



aVx 3 ' 



y=^— (4) 

7 a 2 Vi 312 w 



But since S = d for the conditions under which y is used 

7 = o 2 (5) 

Derivations 

The space charge equation based upon Poisson's equation and the 
energy equation for an electron is 21 

d 2 V 4 

20 In this work e is the charge on the electron in e.s.u., m its mass in grams, c is 
the ratio of electrostatic and electromagnetic units, 1 volt = (10 8 /c) e.s.u., 1 ampere 
= (c/10) e.s.u. 

21 See, for example, Dow, "Fundamentals of Engineering Electronics," John 
Wiley, 1937, pg. 100. 
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This implies that all of the electrons at any plane have the same speed, 
a condition which must be borne in mind when the present analysis is 
applied. 

Integration of equation (6) yields 

S)'-t?<^ + — « (7) 

The choice of zero for this constant leads to type B distributions. A 
negative constant gives type C distributions, and some of the type D 
distributions while a positive value gives type A and type D distribu- 
tions. These will be considered in detail in the sections which follow. 
Transit time solutions are obtained by writing the energy equation 
for an electron in a conservative field which in practical units is 

dx / 2eVl0 8 V' 2 = 3a yin (g) 



dt \ mz ) c 

Solving for / and introducing numerical values 

t = C ± V- 1 ' 2 dx = 1.68 X lO" 8 fv- U2 dx (9) 

= i.68 x io- 8 f (^)~ l y~ U2dv - 

Specialization of this equation for the various types is carried out 
below. 

Integration Constant Zero — Type B 

If the constant in equation (7) is set equal to zero, the next integra- 
tion gives Child's equation, applicable to the type B distribution when 
the correct values of current are used, corresponding to conditions 
before and after the potential zero. Before the potential zero the 
total current, i.e., the arithmetic sum of injected and reflected currents, 
is (2 - Z)I so that 

(2-Z)7=^, . (10) 

where V is the potential x centimeters before the zero. Hence meas- 
uring distance from the first plane in units of S and expressing po- 
tential in units of Vu 

_ __1 ^!_ fin 

a -~(2-zy 2 (2 - zy 12 ' K ' 
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Similar analysis for the region beyond the potential zero where 

a 2 V 3 ' 2 

ZI = -^r> (12) 

yields 

1 <p 3li 

°+ = (2 — ZY 12 Z^* ' 

Introducing the identity y = a 2 in equation (13) gives the relationship 

z v 2 



Zy = 



2 - Z 



<P*i*J. 



+ <P m ■ (14) 



Some of the limiting curves associated with the B solutions are 
closely related, as is implied by the use of a common letter. Curve c 
in Fig. 2 corresponds to maximum values of <p for fixed values of a 
and hence of y = a 2 . By inspection this is seen to be the same condi- 
ion as is implied by curve c in Fig. 9. To find these curves we make 
<p a maximum in (13) with respect to Z holding a constant and obtain: 

M \Z = 2^ 2 /(l + <P m ), (15) 

W I 0- = (1 + 0,1/2)3/22-1/*. (16) 

From these the various other forms of curve c are readily found by the 
relationships y = a 2 and Zy = Za 2 . The curve b of Fig. 2 corre- 
sponds to Z = 1 and gives in equation (13) 

(b) «r = 1 + *,»/«. (17) 

The minimum transmitted current for fixed <p, curve / in Figs. 8 and 9, 
is seen to correspond to y — > oo ; for this the virtual cathode recedes to 
the first plane and 

(0 Zy = ?*. (18) 

Introducing dV/dx obtained from equation (7) with the correct 
current and a zero constant into equation (9) and integrating gives 

cV m 

T = (2 - zyiw 2 + const -' (19) 

which in units of t measured from the first plane with potentials in 
units of Vi gives 

1 - <p'/ 4 
T - " (2 - Z) 1 ' 2 " ( 2 °) 
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Similar analysis to points beyond the potential zero yields 

- 1 -4. £!1 (2i\ 

r + - ( 2 _ zy 12 "*" z 1 ' 2 ' k ; 

Integration Constant Negative — Types C and D 

Introducing for the constant in equation (7) a negative value, say 
— (aVi) 112 , will give a positive value of V (equal to a Vi) for dV/dx = 
with d 2 V/dx 2 > and must therefore lead to solutions of the C type. 

Integrating once more and introducing the unit S gives 

X = =fc So(<p 112 + 2a l/2 ) V 2 - « 1/2 + const. (22) 

Expressing distance from the first plane in units of So, we find two 
possibilities: 

„ D = + ( v i/* + 2a 1 ' 2 ) V 2 - « 1/2 - (1 + 2a 1 ' 2 ) VI - a 1 ' 2 (23) 

and 

a _ = _ (^1/2 + 2a 1 ' 2 )^<p 1 ' 2 - a 1 ' 2 + (1 + 2a 1 ' 2 ) VI - a l/2 . (24) 

The first of these solutions gives a potential distribution rising con- 
tinuously as a increases from zero, hence of type D as was anticipated 
by the subscript. The second solution decreases to a minimum at 

U. = (1 + 2a'/ 2 )Vl - « l/2 (25) 

and then increases, the equation to the right of the minimum being 

a+ = (^,1/2 + 2 «i/2) V 2 - a 1/2 + (1 + 2a 1 ' 2 ) Vl - a l/2 . (26) 

Curves given by equations 23, 24 and 26 are drawn in Fig. 3. If 
values of <r and <p corresponding to conditions on the boundary planes 
are entered in the figure, a C solution is indicated only if the point 
falls upon a curve of the <r+ type. This curve then gives the potential 
distribution to the right of the minimum ; to the left of the minimum 
the distribution is given by the o-_ curve with the same value of a, 
which has the interpretation a = <p min . for this case. Points entered 
on the a- or <x D curves will clearly give D type solutions. 

Three equations for limits of the C region can easily be written 
down on the basis of the above equations : 

(b) c = 1 + <p™. (27) 

(d) <r - (1 + 2^' 2 )V1 - ^ l/2 - (28) 

(e) a = W 2 + 2)V^ 1/2 - 1: (29) 
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The curve a is obtained by making <r + a maximum with respect to a 
while holding <p constant. This gives 

,. f « - «Wn. = *(1 + *> 1/2 )- 2 , (30) 

w U = (1 + ^w. (31) 

Here ^> and <r are coordinates of a point on the a curve and a = v3 rain . 
is the parameter value for the potential distribution curve tangent to 
curve a at that point. The <r + curves give type C solutions for values 
before the tangent point and give C overlap solutions beyond this 
point. 

All the curves described in this section are readily transformed to 
current voltage plots by the relationships 7 = Zy = a 2 . 

The transit times for the various curves are found from equations 
(7) and (9) using the value — (a7i) 1/2 for the constant. Integrating 
and measuring time from the first plane, they are 

TD = + (pi* - a 1 ' 2 ) 1 ' 2 - (1 - a 1 ' 2 ) 1 ' 2 . (32) 

r c _ = - (p 1 ' 2 - a" 2 ) 1 ' 2 + (1 - a 1 ' 2 ) 1 ' 2 . (33) 

t c+ = + (ip 1 ' 2 - a 1 ' 2 ) 1 ' 2 + (1 - a 1 ' 2 ) 1 ' 2 . (34) 

Integration Constant is Positive — Type D 

Type D solutions include those given by equations (23) and (24). 

Other solutions are obtained by giving the integration constant of 
equation (7) a positive value, say + (iSFi) 1 ' 2 . Integrating the equa- 
tion and measuring distances from the first plane in units of S we 
obtain the two possibilities: 

t = + {v m - 2/3"*)V 2 + /3 1 ' 2 - (1 - 2/3 1 ' 2 )Vl + /3 1 ' 2 , (35) 
which applies for ip > 1, and 

a = - W 2 - 2/3 1 ' 2 )Vv 1/2 + /3 1/2 + (1 - 2/3 1 ' 2 )Vl + 0" 2 , (36) 

which applies for ip < 1. Corresponding transit times are: 

r = + (v 1/2 + /3 1/2 ) 1/2 - (1+ W 2 . (37) 

r = - (<p 112 + W* + (1 + /3 ,/2 ) ,/2 . (38) 

Integration Constant is Positive — Type A 

The potential distribution curves of the A type are identical in 
form with those of the D type, where (p < 1, which result from the 
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positive value of the constant in equation (7). They differ in nu- 
merical values in that the current / must be replaced by the value 21 
to allow for the reflected current. The correct equation is then 

a = [(1 - 2/3 1 ' 2 )(l + W - (<p 112 - 2P' 2 )(<p l i* + ^/2)i/2]2- 1 ' 2 . (39) 

The corresponding transit times are given by 

T = [(1 + /3i/2)i/2 - (^1/2 4. ^itf-p-m . (40) 

To the right of <p = the space is free of charge so that the potential 
gradient is constant and equal to the value at <p = obtained by taking 
the derivative of equation (39). This value is 

gg _ 4V2>" (41) 

da 3 

Concerning Completeness 

We may now review our work and see that no possible space charge 
distributions can have been omitted. Starting from the fundamental 
equation (6) we obtain equation (7) with an undetermined integration 
constant. Setting this constant equal to zero we could integrate once 
more, obtaining a solution formally identical with Child's equation. 
If we supposed that the cathode plane defined by the Child's solution 
lay to the right of the initial plane, then the only freedom left in the 
solution was represented by Z, the fraction of current passing through 
the plane. All physically sensible values of Z, i.e., to 1, are included 
in the solutions. If the cathode is assumed to lie to the left of the 
plane, then Z must equal 1 and the solution which arises is given by 
a = in equations (23) or (26), or |8 = in equation (35) — that is, 
a D solution. For a negative value of the constant, a further inte- 
gration gave only two possibilities. Each of these was investigated 
for all possible values of the constant. A similar statement is true 
for positive values of the constant. 

As was stated in the text, space charge distributions corresponding 
to injection from both bounding planes can be handled in formally 
the same way as injection from one plane ; therefore we may conclude 
that all solutions to the problem given by specifying the boundary 
conditions on two planes, subject to the assumptions represented by 
equation (6), have been determined. 
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Equations for Boundary Curves 

For convenience we append a table of equations for the limiting 
curves occurring in the figures. 



Symbol 


Equation 
numbers 


a 


7 


Z 


27 


a 
b 




31 
17.27 

IS, 16 
28 
29 
18 


(1 + ¥>>/»)V» 
1 +<(>»* 

(1 +«»W)W2-W» 


(1 + «>" 1 )' 
(1 + ¥>"«)* 

(1 + v" l )'/2 

(1 +2v>>«)"(l - 

(v>"» + 2)»(*" 1 
oo 


- 1) 


1 
1 

2*>'/» 

(1 +v ,l/2) 

1 
1 




(1 + „ifl)J 
(1 + «*•>' 

V"V(1 + ¥>''»)' 
(1 +2*>/i)'(l - 
( v '/» + 2)«(<p»» 




d 


(1 + 2*>"»)Vl - »>>« 


„ifl) 


t 
1 


(?>»/» + 2W<p" 1 - 1 

co 


- 1) 



The further relationship 

« = -Pmin. = *(1 + ^ 1/2 )" 2 

holds for the a curve. 



(30) 



